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We improve the resummation calculations in the ResBos program for the Higgs boson production 
via gluon-gluon fusion by including the NNLO Wilson coefficient functions and G-functions. The 
improvement leads to increasing the total cross section predictions of the new ResBos program, 
dubbed as ResBos2, for Higgs boson production by about 8% and 6% at the Tevatron and the LHC, 
respectively, as compared to the old ResBos program. Furthermore, the improved predictions are 
compared with those from the programs HNNLO and HqT2. We find that they agree well for the 
CI total cross sections but differ slightly for the transverse momentum Qt distributions. With ResBos2, 

we present the distributions of the two variables ar and <f>* , which can have better experimental 
resolutions than Qt, for the process of Higgs boson decaying into a photon pair. Theoretical 
uncertainties of the ResBos2 predictions are also discussed. 

> 

O ' PACS numbers: 12.38.Cy,12.38.Qk,14.80.Bn 

I. INTRODUCTION 

O |. In the standard model (SM), the Higgs boson is responsible for the origin of clcctroweak symmetry breaking and 
JL ' the generation of elementary particle masses. In the past decades, the Fermilab Tevatron has devoted great effort to 
qj^', searching for Higgs boson in the SM, whose mass is not theoretically precisely predicted. Up to now, the Tevatron 
r~| ■ has excluded a wide range of the Higgs boson mass. With the data corresponding to an integrated luminosity of 10 
fb _1 , the CDF and D0 collaborations found about 2a excess in the mass region 115 GeV < mn < 135 GeV [l[. 

At the CERN Large Hadron Collider (LHC), the ATLAS and CMS collaborations found an about 3a excess 
for the Higgs boson mass in the range 124 GeV < mjj < 126 GeV using the data collected in 2011 at a center-of-mass 
(CM) energy of 7 TeV. Most recently, as a milestone of particle physics, a new resonance with a mass of about 125 
GeV, a Higgs-like particle, has been observed at a 5a level [3, [j| after the data collected at 8 TeV are taken into 
CO , account. In the future, when more data are accumulated, it is promising that LHC can tell whether this is the SM 
\l ■ Higgs boson or something else. 

I/*-) \ The SM Higgs boson is produced at the LHC mainly through the gluon-gluon fusion process. At the lowest order, 
■ the Higgs boson couples with gluons via a top quark loop |6j, and the higher order corrections involve multi-loop 
CNI ' diagrams. The QCD next-to-leading order (NLO) correction for this process has been investigated in both the infinite 
0, H| and finite j^, [To[ top quark mass limits, and can enhance the total cross section by more than 80% for a 125 
GeV Higgs boson produced at the LHC with a CM energy of 7 TeV. The QCD next-to-next-to-leading order (NNLO) 
correction for the inclusive total cross section has been obtained in the infinite top quark mass limit (11H13J, and 
increases the NLO cross section by about 25%. The validity of taking this infinite top quark mass limit at NNLO has 
been studied in Refs. EMI. Furthermore, the soft gluon contributions have been resummed up to next-to-next- 
to-leading logarithm (NNLL) [13], leading to an additional increment of cross section by about 7% at the LHC. The 
next-to-next-to- next-to-leading logarithm (NNNLL) resummation has also been studied in Refs. [l8rl20l |. Recently, 
the contributions in the form of (Carols) 11 , where a s is the strong coupling constant and Ca(= 3) is the Casimir 
operator for the adjoint representation of SU(3) group, have been resummed to all order in the soft-collincar effective 
theory [2ll l22j . Besides the total cross section, a precise prediction of the transverse momentum distribution is of 
crucial imp ortance to understand the properties of the SM Higgs boson, which has been investigated at the QCD 
NLO 1 |23H29| . Moreover, contributions at 0{a b s ) to the production of Higgs boson associated with two jets have also 
been accomplished in Rcf. [30(. Nonetheless, in order to obtain reliable theoretical prediction in the low transverse 
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1 In this paper, for the transverse momentum distribution of Higgs boson production, the QCD NLO results are referred to as those up 
to 0(at). 
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momentum region, the logarithmically-enhanced contributions due to soft gluon emission have been resummed up to 

nnll mm!]. 

In this paper, we show the improvement on the transverse momentum resummation calculations in the RcsBos 
program |43f . in which the large logarithms of the form In" (m 2 H j 'Q^f) have been resummed to all orders in a s . The 
improved ResBos program includes the Wilson coefficient function up to NNLO and the Sudakov exponent with NNLL 
accuracy, and is matched to the C( a s) fixed-order prediction in the large transverse momentum region. Hereafter, 
we refer to the improved version of RcsBos as ResBos2 [44| . We shall compare the RcsBos2 predictions on the total 
cross section and transverse momentum distribution with those from the parton level Monte Carlo programs HNNLO 
H HI] and HqT2 [HM Ej]. The HNNLO program gives prediction exactly at NNLO. The HqT2 prediction 
includes resummation effects up to NNLL order, and is matched to NLO results at large transverse momentum 
region. Though both HqT2 and ResBos2 calculations are performed at the same NNLL+NLO order in predicting 
the transverse momentum distribution of the Higgs boson, the integrated rate of the HqT2 calculations is required 
to reproduce the QCD NNLO total cross section, while there is no such requirement in ResBos2 calculations. The 
difference in their numerical predictions will be discussed in detail in Sec. IIIII 

In the transverse momentum resummation scheme, the values of non-perturbative coefficients, needed for predicting 
the low transverse momentum distribution, arc obtained by fitting to the experiment data of Drcll-Yan processes. 
The non-perturbative coefficients for gg H arc assumed to be scaled by the ratio Ca/Cf from those for Drcll-Yan 
processes, where Cf(= 4/3) is the Casimir operator for the fundamental representation of SU(3) group. However, the 
non-perturbative physics may not follow this assumption, so when we provide the soft gluon resummation prediction 
on gg — > H, the theoretical uncertainty due to the dependence on non-perturbative coefficients is also investigated in 

sec. urn 

This paper is organized as follows. In Sec. [Til wc briefly review the formula of transverse momentum resummation 
for Higgs boson production, and give explicit expressions of NNLO Wilson coefficient functions included in ResBos2. 
In Sec. IIIII wc present the numerical results for the total cross section and the transverse momentum distributions 
given by ResBos2, and compare them with the HNNLO and HqT2 predictions. We conclude in Sec. [TV] 



II. RESUMMATION FORMULA 

In this section we briefly review the Collins-Soper-Sterman (CSS) [45^47| formula of the soft gluon resummation for 
Higgs boson production and decay (40| used in the ResBos program [43 1. By using the narrow width approximation, 
the Higgs boson production process can be factorized from its subsequent decay. The resummation formula for the 
inclusive differential cross section of the gg —> H process can be written as [32|, S3] 

dajgg -» HX) Q 2 Q 2 T H /m H 

dQ 2 dQ\dy Ka ° S (Q 2 - m\Y + (Q 2 T H /m H ) 2 [ > 

{ 1 



(2tt)s 



d 2 be iQTb W gg (b*,Q, Xl ,x 2 ,C lt2 , 3 )W" p (b, Q, x u x 2 ) + Y{Q Tl Q, x u x 2 ,C A ) }, 



where S is the square of the CM energy; Q, Qt, y, and Tjj are the invariant mass, transverse momentum, rapidity, 
and total decay width of the Higgs boson, respectively. In Eq. ([1]), the quantity <tq, defined as <t lo (oo, 0, 0), arises as 
an overall factor in the infinite top quark mass limit. It is given by 
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V2G F a 

1 576tt ' v ; 

where Gf is the Fermi constant and a s is evaluated at the hard scale C 2 Q. Furthermore, we have multiplied tro by 
an additional factor 

_ <7 LO (m t ,m b ,m c ) 
K ~ a L O(oo,0,0) ' [) 

in order to take into account the mass effects of the top, bottom, and charm quarks at LO. It has been shown that 
multiplying the NLO cross section in the infinite top quark mass limit by the factor k is a good approximation to the 



full mass-dependent NLO cross section over a wide range of the Higgs boson mass [48( . The total cross section can 



be obtained by integrating the result in Eq. (fTJ) over the transverse momentum and the rapidity. 
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In Eq. ([T]), the term containing W gg dominates at small Qt, and behaves as Q T 2 times a series of In™ (Q 2 /Q 2 ^) 



The explicit form of W gg can be expressed as 



W gg (b,Q,x 1 ,x 2 ,C 1 , 2t3 ) =e - s ^ c ^ {C 3 a®fa)(xi)(C gb ®f b )(x 2 ), 



(4) 



where the Sudakov exponent is given by 



S(b,Q,d,C 2 ) = 



C'lQ 2 



ClQ 2 



+ S(a s (/2),d,C 2 ) 



^(a s (/2),d)ln 

The coefficients A and -B and the Wilson coefficient functions C ga can be expanded as a power series in a s : 

'MAT 



A{a s (p),C x ) = 

n=l 
oo 

B(a s (p),Ci,C 2 ) = 



7T 

a s {p.) 



B (n) (Ci,C 2 ) 



(5) 

(0) 
(7) 



and 



C ga (z,b,fi,C 1 ,C 2 ) 



(8) 



Here, the constants Ci and C2 are included in order to adjust the scales at which the coefficients A and B are 
evaluated, so that the higher-order corrections are moderately small [471 ]. The renormalization scale /1 in C ga is, 
in principle, arbitrary and set equal to C3/6 in the RcsBos program. The constant C4 in the non-singular part 
Y of Eq. ([T]) sets the scale of the fixed-order perturbative calculation which is of the order of the Higgs boson 
invariant mass Q. Unless specified otherwise, we use the canonical choices for the renormalization constants, i.e., 
C\ = C3 = bo = 2e~ lB ,C 2 = C4 = 1, in our numerical calculations. The perturbation expansion quantities in Eqs. 
©-© can be extracted order-by-order from fixed-order calculations. In the old version of the ResBos program, we 
have included A^' 2 ' 3 \ £?( 1,2 ) and C^ 0,1 -', whose analytical expressions are given in Ref. [ioj . 

In this paper, we improve the ResBos calculations by including the NNLO Wilson coefficient functions, which can 
be readily derived from 'H. g ^Sz ab , as given in Eqs. (23) and (24) of Ref. [49|. As shown below, the explicit expressions 
of the NNLO Wilson coefficient functions are 
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Here rif is the number of active quark flavors, and £ 3 is the Riemann zeta function which is defined as 
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The above coefficients can be derived from Eq. (33) of Ref. [49|, 
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(12) 



where <J a b^H is the partonic cross section of Higgs boson production process ab — > H, and Qq is a small cut-off scale 
of transverse momentum, which is taken to be much less than m# . The left hand side of Eq. (|12[) can be expressed 
as 



d 1r —^T^-(lT,rn H ,s = m 2 H /z) = z er Rab^H(z,m H /Q ) , 



(13) 
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in which Co is the Born level partonic cross section given in Eq. ([2]), and 

*3L*(*, m H /Q ) = it eS52w + n e» w + ii =222 w + « (=^2 w - i6o,x£S2 w) 

+ (wf^W - 4Cs Ef^C*)) + ©(Qg/mir) , (14) 



with Iq = log(m^/Qo) S3- ^fg^-2b^ n = •"•) •••>^ are known in the transverse momentum resummation scheme 37|. 
In the limit Qo — > 0, the terms 0(Q^/m 2 H ) can be neglected. Moreover, 

Wf^o* (*) = ^C-^(z) + <^C$(z) + C« ® c£>(z) + GW ® G${z). (15) 

The calculation of W^_ (,(z) is challenging and has been originally performed in Rcf. [49j , We have numerically 
checked Eq. (|12[) independently using the analytical expressions for the total cross section (the first term on the right 
hand side of Eq. (|12p) in Ref. fl2j and transverse momentum distribution (the second term on the right hand side 
of Eq. (|12p) in Ref. (2(|, and found some typos for H^g<^ gq in the first version of Ref. [49| 2 , while the expression of 
Uggp-gq in the program HNNLO [H, H| is correct. 

The function W^L P describes the non-perturbative part of the soft gluon resummation. We use the BLNY param- 
eterization form [5(|, with the non-perturbative coefficients scaled by the factor of Ca/Cf = 9/4. This scaling factor 
is applied since the initial states are gluons in the gg —> H process, instead of quarks in the Drell-Yan process. The 
effect of varying the non-perturbative coefficients will be discussed in the next section. 

In addition, the G-function, which is only present for the gluon initiated processes, can be incorporated into the 
CSS formula as shown in Ref. [4!|. Therefore, the term W can be further improved as 

W gg (b,Q,x 1 ,x 2 ,C lt 2, 3 ) = e - s ^ Q > c ^ Y, [( C 9* ® fa)(xi)(C gb ® f b )(x 2 ) + (G ga <g> / )(x 1 )(G a6 ® fb){x 2 )] • 

a,b=q,q,g 

(16) 

The lowest order of G-function is at 0(a s /7r), and 

oo 

G ga (z,a s ) = -G«(z)+^p) GH(z), (17) 

where the coefficient functions are known as [4!| [5l|, [52[ 

G V{z) = C A ^, G W(z) = C f ^, (18) 

in which z denotes the fraction of momentum carried away by the gluon. 

Finally, in Eq. (fl[). the term containing Y incorporates the remainder of the cross section which is not singular 
as Qt — > 0. It consists of the difference between the full fixed-order cross section and the small-Qr limit of this 
cross section. At small Qt, they cancel with each other, so that the contribution of the Y-term is small, and the W 
term dominates. At large Qt, where the logarithmic terms become small, the W term cancels against the small-Qr 
limit term (to the given order in a s ), so that the result approaches the fixed-order calculations. More details of how 
this matching process between the resummation calculation and the fixed-order one is implemented in the ResBos 
program can be found in Ref. (43j . 

III. NUMERICAL RESULTS 

In this section, we investigate the effect of including the NNLO Wilson coefficient functions and G-function con- 
tributions in the RcsBos2 calculations on the total cross sections and transverse momentum distributions. Unless 
otherwise specified, the parton distribution function (PDF) used in the numerical calculations is chosen as CTEQ6.6 



2 After our paper appeared on the arXiv, the expression of H^n^lqq in Ref. [49^ was corrected in its revised version 
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|53(. The scale uncertainty is obtained by varying the scale C2 around the canonical choice (C2 = 1) by a factor of 
two, but maintaining the relations 

d = C 2 b , C 3 = b and d = C 2 . (19) 

Keeping these relations between C,,i = 1,2,3,4, the Wilson coefficient functions are not changed. Therefore, the 
dominant effect of the variation is to change the shape, rather than the overall size of the cross section. We shall 
compare the ResBos2 predictions with those from the programs HNNLO and HqT2 on the total cross sections and 
the transverse momentum distributions. 





m„ (GeV) 


ResBos2 


ResBos 


HNNLO (NNLO) 


HqT2 (NNLL+NLO) 


Tevatron 


115 




U - J1 -6.9% 


0.96lf 3 ; 6 7 * 


g7 +14.2% 


120 


0-87™ 


n on+15-8% 
U - 8U -6.6% 


0.84™ 


n or + 13.9% 
U - 00 -13.8% 


125 


77+ 9 - 1% 
'-6.5% 


q 71 + 15.9% 




n 7 c + 13.8% 
U -' J -13.8% 


130 


n fi o+10.3% 
U - b8 -5.9% 


n fi o + 15.9% 
u - Di5 -6.4% 


0.66+^1 




LHC 7 TeV 


115 


1 c 1 1 +7.8% 
10 - 1J --5.8% 


1 4 91 + 8 - 3% 


15 l6 +9 - 0% 
io.io_ 10 7% 




120 


1 Q 00 + 7.8% 

i0,oy -5.7% 


13.06™ 


1 ■? sn+iO-6% 
13.»U_ 10 5% 


13 94+10-2% 
io.yi_ 10 4% 


125 




12.03+*;^ 


1 9 79+10.2% 

iz -' z -10.6% 


1 9 CQ+10.7% 
iZ -°°-10.4% 


130 




n -, 2 +8.6% 


, , 7 r + 10.8% 


11 S4+ 9 - 9% 
J-J--»4_ 10 4% 


LHC 8 TeV 


115 


19 15 +7.5% 

i».±o_ 6 B% 


17.58+^* 


19.0V 


1Q OK + 10.8% 
u.z,o_ gs% 


120 


ll - 00 -6.S% 


16.21+^* 


17 CQ+9.7% 
' — 9.9% 


17 76 +9 - 8% 

1 '-'°_10.1% 


125 


16 3l+ 7 - 5% 

iD - 31 -6.4% 


14.98+^* 


16 26 +10 2% 
ro.zo 10 4% 


16.39 +9 8 ^ 


130 


15 11 +7 - 5% 

i0 - 11 -6.4% 


13.89™ 


1 c n7 +10.9% 
i0 - U ' -10.6% 




LHC 14 TeV 


115 


48.84+ 7 ; 6 7 | 


46 n3+ 10 - 0% 


49 24+ 9 ' 4% 

^y-^i_9.8% 


48.901 9 ; 9 ! 


120 




42.84± 9 ' 8 ^ 


45 57 +io.4% 

* -°'-9.5% 


45 52 +10 ' 3% 

40-0^-8.4% 


125 


42.39l£S* 


39 96 +9 ' 6% 


42 61 +9 6% 


42 57+ 9 - 7% 

4z -0'_8.8% 


130 


39.651™* 


17 QS+9- 4 % 
J '- JO _4.8% 


oq qq+ii-1% 
oa.yo_ g g% 


09 75 +9.7% 



TABLE I: The ResBos2, ResBos, HNNLO and HqT2 predictions on the total cross sections (in pb) for Higgs boson production 
via g + g -> H + X at the Tevatron (1.96 TeV) and LHC (7 TeV, 8 TeV and 14 TeV). The upper (lower) uncertainties, expressed 
in the form of percentages, are obtained by dividing (multiplying) the canonical scale by a factor of two. 

Table U shows the predictions on the total cross sections for Higgs boson production at the Tevatron and the LHC. 
By including the NNLO Wilson coefficient function C^ 2 \ ResBos2 predictions at the Tevatron and the 7 TeV (8 
TeV, 14 TeV) LHC are larger than ResBos predictions by about 8% and 6% (9%, 6%), respectively. The theoretical 
uncertainties of ResBos (RcsBos2) predictions are estimated by varying the scales in the resummation formula, as 
discussed around Eq. (fTTJ)) . They are reduced in RcsBos2 to about 9% and 8% (7%, 8%) at the Tevatron and the 7 
TeV (8 TeV, 14 TeV) LHC, respectively 

We can also see that RcsBos2 predictions agree with the predictions from HNNLO and HqT2 programs within a 
couple of percent. When the CM energy of the LHC is increased from 7 TeV to 8 TeV, the total cross sections for 
the Higgs boson production are increased by about 27% for m#=125 GeV. 





CTEQ6.6 


CT10 NLO 


CT10W NLO 


CT10 NNLO 


MSTW2008NNLO 


NNPDF2.3NNLO 


Tevatron 


0.77 ± 7.2% 


0.77 ± 7.2% 


0.76 ±7.1% 


0.77 ± 7.2% 


0.78 ± 6.6% 


0.80 ± 4.5% 


LHC 7 TeV 


12.80 ± 3.3% 


13.33 ± 3.3% 


12.82 ± 3.3% 


12.65 ± 2.9% 


12.69 ± 2.7% 


13.73 ± 2.3% 


LHC 8 TeV 


16.31 ± 3.4% 


16.53 ± 3.4% 


16.95 ± 3.4% 


16.63 ±3.1% 


16.30 ± 2.4% 


16.90 ±2.1% 


LHC 14 TeV 


42.39 ± 3.5% 


42.64 ± 3.5% 


42.91 ± 3.4% 


41.87 ±3.6% 


43.10 ±2.3% 


43.28 ± 2.0% 



TABLE II: The total cross sections (in pb) for Higgs boson production via g ± g — > H ± X at the Tevatron (1.96 TeV) and 
LHC (7 TeV, 8 TeV and 14 TeV) by using different PDF sets in ResBos2. The PDF induced uncertainties are estimated at (or 
rescaled to) 90% confidence-level, and expressed in the form of percentages. 

To study the dependence of ResBos2 predictions on the choice of PDF sets, we show in Table |TT] the total cross 
sections of producing a 125 GeV Higgs boson at the Tevatron and the LHC, and corresponding PDF uncertainties when 
using the CTEQ6.6, CT10 NLO, CT10W NLO, CT10 NNLO HHH, MSTW2008NNLO [Hg and NNPDF2.3NNLO 



7 



57 PDF sets. We compare these widely used PDF sets as recommended by the PDF4LHC Working Group in Ref. 
58|. There is an about 8% difference between the NNPDF2.3NNLO and CT10 NNLO predictions at the 7 TeV LHC. 



In general, the NNPDF2.3NNLO sets predict larger cross sections for the Higgs boson production than the other PDF 
sets. A similar conclusion is also found in Ref. 1591. 
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FIG. 1: The different theoretical predictions on the transverse momentum distributions for the Higgs boson production at the 
Tevatron (1.96 TeV) and the LHC (7 TeV, 8 TeV and 14 TeV). In the bottom of each plot, the ratios to ResBos2 predictions 
are also shown. 



In Fig. [I] we present ResBos2 predictions on the transverse momentum distributions of the Higgs boson, compared 
with the predictions of ResBos and HqT2. The predictions of RcsBos2 generally increase the overall size of the 
distributions, compared to ResBos. Especially, the improvement in the small transverse momentum region can be as 
large as 20% at both the Tevatron and the LHC. The peak positions are slightly shifted to small Qt region. Meanwhile, 
ResBos2 predictions are smaller than ResBos by about 20% for Qt ~ 60 GeV (80 GeV) at the Tevatron (LHC). This 
behavior is due to the fact that the inclusion of NNLO Wilson coefficient functions modifies the resummed (W) 
contribution and hence the matching condition in the intermediate Qt region. We note that the high Qt predictions 
remain to be the same. The HqT2 predictions at the LHC (7 TeV, 8 TeV and 14 TeV) in small Qt region are 
almost the same as ResBos2, but the peak height of the HqT2 prediction at the Tevatron is a little lower than that 
of RcsBos2. On the other hand, the HqT2 predictions at moderate Qt region are closer to the ResBos predictions, 
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and can be about 20% higher than RcsBos2 predictions for Qt ~ 60 GeV (80 GeV) at the Tevatron (LHC). 
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FIG. 2: The transverse momentum distributions of the Higgs boson for different scale choices at the Tevatron (1.96 TeV) and 
the LHC (7 TeV, 8 TeV and 14 TeV). 



Figure [2] shows the resummation scale dependencies of the transverse momentum distributions of the Higgs boson 
produced at the Tevatron and the LHC. To estimate the resummation scale dependencies, we vary the hard scale 
coefficient Ci by a factor of two around the canonical choice, but hold the relations between (7j, i = 1, 2, 3, 4, as shown 
in Eq. ([T9"]) . With this choice, the Wilson coefficient functions C gg and C gq do not change when Ci varies. Therefore, 
mainly the shapes of distributions are affected by varying the resummation scales. As shown in Fig. [2] the shapes 
of the transverse momentum distributions of Higgs boson are changed at both the Tevatron and the LHC. The peak 
positions of the distributions are shifted by several GeV's when the scales are varied. Generally, the distributions in 
the low and high Qt regions are increased and suppressed, respectively, when C*2 increases. 

Figure [3] shows the ratios between the transverse momentum distributions with and without G-functions at the 
Tevatron and the LHC. It can be seen that the effect of G-functions is very small. As shown in the figure, by including 
G-functions, the shapes of transverse momentum distributions are slightly changed, by less than 1%. This small effect 
is understood, for the contributions of G-functions to the transverse momentum distribution of the Higgs boson start 
at 0(a 4 s ) [sec Eqs. ((TB I) - (fig ]) ]. Notc that thc L0 cross section is at 0{a 2 s ). 

Figure |4] presents the dependencies of the transverse momentum distributions on non-pcrturbative coefficients in 
soft gluon resummation at the Tevatron and the LHC. Based on the best-fit non-perturbative coefficients in the 
BLNY parameterization form [H(J, but scaled by the ratio Ca/Cf = 9/4, the dependencies on thc non-perturbative 
coefficients are shown by multiplying or dividing the best-fit values by a factor of two. We can see that when the 
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FIG. 3: The ratios between the transverse momentum distributions with and without G-functions at the Tevatron (1.96 TeV) 
and the LHC (7 TeV, 8 TeV and 14 TeV). The oscillations of the ratio curves in the figure are due to numerical uncertainties. 



values of the non-perturbative coefficients increase, the peak positions are slightly shifted towards larger Qt at both 
the Tevatron and the LHC, and the peak heights are decreased by a few percent at the Tevatron but hardly change 
at the LHC. On the other hand, when the values of the non-perturbative coefficients decrease, the peak positions 
slightly move towards smaller Qt, but the peak heights remain the same at both the Tevatron and the LHC. Hence, 
we conclude that the transverse momentum distributions of the Higgs boson produced at the LHC are dominated 
by the perturbative Sudakov resummation effect of multiple soft gluon emissions, and the effects of non-perturbative 
physics are not as important. 

In the very low Qt region, the non-perturbative physics becomes dominant, whose effect in the prediction of the 
Higgs boson transverse momentum distributions can be better estimated after the non-perturbative parameters (and 
the function form) are measured from the gluon initiated scattering processes, such as top quark pair production at 
the LHC. 

Before closing this section, we would like to discuss the prediction of ResBos2 program on the kinematical distribu- 
tions of the two photons from the decay of Higgs boson produced at the Tevatron and the LHC. When measuring the 
transverse momentum distribution of the Z boson produced at the Tevatron, the experimentalists have proposed two 



other variables, i.e. cit [6fJ and <fi* [61(, which have been shown to be less susceptible to the effects of experimental 



resolutions and efficiencies. For example, the D0 collaboration at the Tevatron has measured the distribution of 
4>* for the Z boson production and its decay into two leptons [f32j]. In the process ,917 — > H — > 77, it is convenient 
to measure these two variables, similar to what was done in the analysis of the Drell-Yan lepton pair production. 
Without restating their definitions given in Rcfs. [r30l. l6lj. we show the normalized distributions of 4>* and ar for 
different transverse momentum cuts on the individual photon in Fig. O and Fig. [5J respectively. It is evident that 
the distributions of 4>* are almost insensitive to the transverse momentum cuts on the individual photon at both the 
Tevatron and the LHC. From Fig. [6j we can see that the distributions of clt change a little when the transverse 
momentum cut on the individual photon varies from GeV to 50 GeV. Figure [7] and Figure [5] show the scale de- 
pendencies of the distributions of (jf and or, respectively. It can be seen that in the small value region, the scale 
uncertainties of the distributions of 4>* and ax is about ±5% and ±8% at most, respectively. In the large value region, 
the scale uncertainties are negligibly small. 



IV. CONCLUSION 



We have improved the resummation calculations in the ResBos program for the Higgs boson production via gluon- 
gluon fusion by including the NNLO Wilson coefficient functions and G-functions. We show that including NNLO 
Wilson coefficient functions increases the total cross section predictions of ResBos for a 125 GeV Higgs Boson pro- 
duction by about 8% and 6% at the Tevatron and the LHC, respectively. The predictions of the newly improved 
ResBos progran, dubbed RcsBos2, agree with the predictions of HNNLO and HqT2 programs within a couple of 
percent at the Tevatron and the LHC. We also show that for the transverse momentum distributions, the predictions 
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FIG. 4: The transverse momentum distributions for different non-perturbative parameter choices at the Tevatron (1.96 TeV) 
and the LHC (7 TeV, 8 TeV and 14 TeV). 



of ResBos2 generally increase the overall size of the distributions in the low to medium transverse momentum regions, 
as compare to the old ResBos program. Especially, the improvement in the small transverse momentum region can 
be as large as 20% at both the Tevatron and the LHC, and the peak positions are slightly shifted toward lower values 
than ResBos. At the LHC, RcsBos2 predictions for the transverse momentum distribution are almost the same as the 
HqT2 program. At the Tevatron, the peak height of RcsBos2 prediction is somewhat higher than that of the HqT2 
program. When increasing the renormalization constant C2 and keeping the relations between Cj, i = 1, 2, 3, 4, in the 
ResBos2 calculations, the distributions in the low and high Qt regions are increased and suppressed, respectively. The 
inclusion of G-functions in RcsBos2 can mainly slightly modify the shape of the transverse momentum distribution 
of the Higgs boson at the LHC, by less than 1%. Similarly, varying the non-perturbative coefficients by a factor of 
two only slightly modifies the peak position and the height of the Qt distribution. Thus, the transverse momentum 
distributions of the Higgs boson produced at the LHC are dominated by the perturbative Sudakov resummation effect 
of multiple soft gluon emissions, and the effects of non-perturbative physics are not as important. 

From an experimental viewpoint, we also present the distributions of two variables ax and <jf in the process 
gg — > H — > 77, which have better experimental resolutions than Qt itself. We find that the distributions of <fi* are 
almost insensitive to the transverse momentum cuts on the individual photon at both the Tevatron and the LHC. 
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FIG. 5: The normalized (j>* distributions for different transverse momentum cuts (PJ, > 0, 25, 50 GeV) on the individual photon 
at the Tevatron (1.96 TeV) and the LHC (7 TeV, 8 TeV and 14 TeV). 

The scale uncertainties of the distributions of </>* in the small value region are about ±5% at most, and the scale 
uncertainties in the large value region are negligible. Similar conclusions also hold for the ar distributions. 
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